Quantum computation of the Anderson transition in presence of imperfections 
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We propose a quantum algorithm for simulation of the Anderson transition in disordered lattices 
and study numerically its sensitivity to static imperfections in a quantum computer. In the vicinity 
of the critical point the algorithm gives a quadratic speedup in computation of diffusion rate and 
localization length, comparing to the known classical algorithms. We show that the Anderson 
transition can be detected on quantum computers with 7—10 qubits. 

PACS numbers: 03.67.Lx, 24.10,Cn, 72.15.Rn 



The problem of metal-insulator transition of noninter- 
acting electrons in a disordered potential was pioneered 
by Anderson in 1958 [l|. Since then it continues to at- 
tract an active interest of researchers all over the world 
(see e.g. 0, 0, 01 an d Refs. therein). In addition to 
analytical and experimental studies of the problem an 
important contribution to the understanding of its prop- 
erties was made with the help of numerical simulations 
based on various computational methods adapted to the 
physics of this phenomenon. Indeed, the numerical stud- 
ies allowed to obtain some values of critical exponents in 
the vicinity of the transition and to study certain system 
characteristics at the critical point including level spac- 
ing statistics and conductance fluctuations for the cases 
of different symmetries and system dimensions (see e.g. 
0, H 0, 0). These numerical simulations are per- 
formed with the help of modern supercomputers and are 
at the border of their computational capacity. 

The recent progress in quantum computation demon- 
strated that due to quantum parallelism certain tasks 
can be performed much faster on a quantum computer 
(see 0] and Refs. therein). The most known example 
is the Shor algorithm for factorization of large integers 
0] which is exponentially faster than any known classi- 
cal algorithm. A number of efficient quantum algorithms 
was also proposed for simulation of quantum evolution of 
certain Hamiltonians including many-body quantum sys- 
tems llflllll and problems of quantum chaos [12lll3Lll4l |. 
In Ref. [13j it was shown that the evolution propagator 
in a regime of dynamical or Anderson localization can 
be simulated efficiently on a quantum computer. How- 
ever, the algorithm proposed there requires a significant 
number of redundant qubits and is not accessible for an 
experimental implementation with a first generation of 
quantum computers composed of 5 - 10 qubits. 

In this paper we propose a quantum algorithm for a 
quantum dynamics in the regime of Anderson localiza- 
tion. This algorithm requires no redundant qubits thus 
using the available n q qubits in an optimal way. The 
propagation on a unit time step is performed in O(n^) 
quantum elementary gates while any known classical al- 
gorithm requires 0(2 nq ) operations for a vector of size 



N = 2 rlq . Due to these properties the Anderson transi- 
tion can be already detected on a quantum computer 
with 7-10 qubits. The basic elements of the algo- 
rithm involve one qubit rotations, controlled phase shift 
C{4>), controllcd-NOT gate Cm and the Quantum Fourier 
Transform (QFT) All these quantum operations have 
been already realized for 3-7 qubits in the NMR-based 
quantum computations reported in Refs. 

[HO]!. Thus 

the main obstacle for experimental detection of the An- 
derson transition in quantum computations is related to 
the effects of external decoherence [Tj| and residual static 
imperfections |l8j | which restrict the number of available 
quantum gates. The results obtained for operating quan- 
tum algorithms show that the effects of static im- 
perfections affect the accuracy of quantum computation 
in a stronger way comparing to the case of random noisy 
gate errors. Due to that in this paper we concentrate our 
studies on the case of static imperfections investigating 
their impact on the system properties in the vicinity of 
the Anderson transition. 

To study the effects of static imperfections in quantum 
computations of the Anderson transition we choose the 
generalized kicked rotator model described by the unitary 
evolution of the wave function tp: 



jP = Uip = exp(-iV(6, t)) exp(-iH (n))ip 



(1) 



Here tp is the new value of ip after one map iteration 
given by the unitary operator U, Hq(ji) gives the rota- 
tional phases in the basis of momentum h — —id/d6, 
the kick potential V(6,t) depends on the rotator phase 
9 and time t measured in number of kicks, ip(9 + 2ir) = 
ip(6). For V(8,t) = kcos0 and H Q = Tn 2 /2 one has 
the kicked rotator model described in detail in [20| . 
The evolution given by (IJ results from the Hamilto- 
nian H = Ho(n) + V(6,t)Si(t), where 5i(t) is a peri- 
odic <5-function with period 1 and [n, 9) are conjugated 
variables. In the case when the potential V(6,t) = 
— 2 tan _1 (2fc(cos 9 + coscoit + cosuiit)) depends quasi- 
pcriodically on time t the model can be exactl y re duced 
to the three-dimensional (3D) Lloyd model 21] . In- 
deed, the time dependence of V(9, t) can be eliminated 
by introduction of extended phase-space with a replace- 
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FIG. 1: (color online) The time evolution of the probability 
distribution \ip n \ 2 in the localized (left column, k — 1.2) and 
delocalized (right column, k = 2.4) phases for n q — 7 qubits 
(N = 2 n "), with < t < 400 (vertical axis) and -N/2 <n< 
N/2 (horizontal axis); k c = 1.8. The color is proportional to 
probability: blue/black for zero and red/ white for maximal 
values. The strength of static imperfections is e = fi — for 
top row and e = /j, = 10~ 4 for bottom row. 



ment Hq — > Ho(n) + u>ini + ^2^2- Then the linear 
dependence on quantum numbers 711,2 gives fixed fre- 
quency rotations of the conjugated phases 9\^. — oJi^t. 
The extensive studies performed in |2lJ showed that this 
model displays the Anderson metal-insulator transition 
at k = k c w 0.5 with the critical exponents being 
close to the values found in other 3D solid state mod- 
els. In this paper following [2^ we choose in Q the 
potential V(9,t) = fc(l + 0.75 cos w\t cos w^t) cos 9 with 
uj\ = 27rA _1 , lo 2 = 2ir\~ 2 and A = 1.3247... being the 
real root of the cubic equation x 3 — x — 1 = 0. The 
rotation phases Ho(n) are randomly distributed in the 
interval (0, 27r). This model shows the Anderson transi- 
tion at k c w 1.8 with the characteristics similar to 
those of the Lloyd model studied in j2]| . 

The quantum algorithm simulating the time evo- 
lution of this model is constructed in the following 
way. The quantum states n = 0, ...,N — 1 are rep- 
resented by one quantum register with n q qubits so 
that N = 2 nq . The initial state with all probability at 
no = corresponds to the state |00...0 > (momentum 
n changes on a circle with N levels). The phase 
rotation Ut — exp(—iHo(n)) in the momentum basis n 
is performed with the help of quantum random phase 
generator built from two unitary operators U T ^ and 
U T 2) . The operator U T 1] = fT?=i gives rotation 

of qubit j by a random phase cf>j. Here and below 



a x ,<7 y ,a z are Pauli matrices. To improve the indepen- 
dence of quantum phases we then apply the operator 

u t^ = rifeli C N {iM-k,3M-k) Ilfeli e^'^^CNiikJk)- 
This transformation represents a random sequence with 

M one-qubit phase shifts e*^* CT ^* and controlled-NOT 
gates CN(ik,jk) followed by the inversed sequence of 
controlled-NOT gates C N (i M -k, jM-k)- Here C N (i k ,jk) 
inverts the qubit j\ if the qubit ik is 1; ik,jk and 
phases 0'- fc are chosen randomly. The resulting random 

quantum phase generator Ut — U^U^~ gives more and 
more independent random phases with the increase of 
M. We use M « 2n q (at n q w 10) that according to our 
tests generates good random phase values. This step 
involves 3M + n q quantum gates. After that the kick 
operator Uj. = exp(— ik(t) cos 9) is performed as follows. 
First, with the help of the QFT the wave function is 
transformed from momentum n to phase 9 representa- 
tion in 0(n q /2) gates. Then 9 can be written in the 
binary representation as 9/2ir = 0.aia2-.a„ q with dj = 
or 1. It's convenient to use the notation 6 = nai + 9 
to single out the most significant qubit. Then due to 
the relation cos# = (— l) ai cos0 = a\ cos9 the kick 
operator takes the form U k = e -' lk ^ cose = e -^f*(*)«»e ) 
where a[ z ' x ^ act on the first qubit. This operator can be 
approximated to an arbitrary precision by a sequence 
of one-qubit gates applied to the first qubit and the 
diagonal operators S m — e vmai6 . The S— operators 
are given by the product of n q — 1 two-qubit gates 
as S m — n™=2 Ci iJ -(7rm2 - - J+:L ) where controlled phase 
shift gate Cj 1 j 2 ((f>) makes a phase shift e 1 ^ if both 
qubits j'1.2 are 1. Then we introduce the unitary opera- 
tor Hy(6) = HS 1 H e- i- i< HS- 2 H e -^ a * HS 1 H 
where H = (af + crf)/v2 is the Hadamard 
gate. It can be exactly reduced to the form 
i? 7 (£) = cos 2 I - sin 2 |cos(20) - ia\ sin7cos(0) + 
iaf sin ^ s'm(29) and hence for small 7 we have 
Ry(8) = e -^7cose + i a f£ s i n ^2e) + 0(7 3 ). 
The term with 7 2 can be eliminated using the 
symmetric representation R-y/2(@)R-y/2(~ @) 
IIS' II e~ i: k°l HS- 2 H e'^"' HS 2 H e^i"" HS^H = 
e -K7cos(e) _^ 0( 7 3). Thus the kick operator is given by 
U k = (R 7/2 (9)R l/2 l-9)) 1 + 0(lj 3 ) where the number of 
steps I — k/j and we used in our numerical simulations 
the small parameter 7 = k/l w 0.2 that gives I w 5 — 10 
for k ~ 1 — 2. After that the state is transfered to 
the momentum representation by the QFT. Thus an 
iteration is performed for 2 n i states in n g elementary 
gates where n g = 2[k/^](n q + 2) + n 2 + 6n q + 3M + 9 
with the square brackets denoting the integer part. This 
algorithm is optimal for the kicked rotator model with 
moderate values of k where n g value remains reasonable. 
It can be easily generalized to d dimensions. 

In our numerical simulations we study the effects of 
static quantum computer imperfections considered in 
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FIG. 2: Top row: logarithm of probability log 10 \ipn\ 2 vs. mo- 
mentum n after t = 10000 iterations; dark gray curves are 
shifted down by 5 (left) and 2 (right). Bottom row: depen- 
dence of IPR £ on time t. The left/right column corresponds 
to localized/delocalized phase at k — 1.2 and k = 2.4 respec- 
tively. The three curves represent e = 0;2 x 10~ 5 ;6 x 10" 5 
with color changing from light gray to black with increase of 
e; /i = e, n q = 10. 

[Til Hil llij . In this case all gates are perfect but 
between gates ip accumulates a phase factor e 1 ^ with 

= Sj(% cr I + Mjffof+i)- Here ^jj/^j var y randomly 
with j = 1, ...,n q , rjj represents static one-qubit energy 
shifts, — e/2 < r?j < e/2, and /Ltj represents static inter- 
qubit couplings on a circular chain, —fi/2 < < /x/2. 

An example of time evolution of probability distribu- 
tion in the momentum representation n is shown in Fig. 1. 
Below the Anderson transition (fc < k c ) the probability 
remains bounded near initial value uq, while above it 
(k > k c ) a diffusive spreading in n takes place. Com- 
paring to the ideal quantum computation the static im- 
perfections lead to probability transfer on levels located 
far away from the center of the wave packet. This ef- 
fect is related to the structure of the QFT where a mis- 
match in the quantum gates generates high harmonics. 
As a result static imperfections create a plateau in the 
probability distribution which level grows with the in- 
crease of e and \i (see Fig. 2). This leads to an artifi- 
cial diffusion of the second moment of the distribution 
< n 2 >=< ip n \(n — tt-o) 2 l^n >• Since the plateau in 
probability extends over all N levels the rate of this dif- 
fusion grows exponentially with n q at fixed e, ufdata not 
shown). A similar effect was discussed in |23j for the 
quantum computation of the kicked rotator with noisy 
gates. Due to that the most appropriate characteristic to 
study is the Inverse Participation Ratio (IPR) £ which is 
extensively used in systems with localization 0, 0] and 
which determines the number of levels on which the wave 
function is concentrated (%2 n |VVi| 4 = I/O- ln contrast 
to < n 2 >, the IPR £ remains stable with respect to noise 
in the gates during polynomially large times [2^ • 

The variation of £ with time and e, ji is shown in Fig. 2. 




FIG. 3: Dependence of the IPR £ and the excitation prob- 
ability W (full and dashed curves for left and right scales 
respectively) on the kick strength k for n q — 10 and t > 10°, 
e = 0;10~ 5 ;2 x 10~ 5 ;4 x 10" 5 ;8 x 10 -5 (corresponding to 
curves from right to left); fi — 0. 

For moderate imperfections, during a rather long time in- 
terval £ remains close to its value in the exact algorithm. 
However, at very large times t > 10 5 it saturates at some 
value which depends on k and e, /i. A typical example of 
such a dependence is presented in Fig. 3. Here, £ shows 
a sharp jump from small (£ ~ 1) to large (£ ~ N) values 
which takes place in a narrow interval of k values. This 
is a manifestation of the Anderson transition from local- 
ized to delocalized states. The critical point k c can be 
numerically defined as a such value of k at which £ is at 
the middle between its two limiting values. The data of 
Fig. 3 show that the critical point k c (e) decreases with the 
increase of the strength of imperfections. The physical 
origin of this effect is related to the additional transitions 
induced by static imperfections which naturally lead to 
a derealization at a lower value of k compared to the 
ideal computation. Another method to detect the posi- 
tion of the critical point k c (e) in presence of imperfections 
is to measure the two most significant qubits which code 
the value of momentum n. After a few tens of measure- 
ments of first 2 qubits one determines the probability 
w = ]C„=(A74,3A74) \i>n\ 2 - At sufficiently large t this 
probability shows a sharp jump from a value W = to 
W w 0.5 when k is varied. This allows to determine the 
critical point and gives the values of fc c (e) close to those 
obtained via IPR £ (see Fig. 3). 

The shift of the critical point Afc c (e) = k c — k c (e) de- 
pends on e, ji and n q . From the IPR data obtained for 
various e, n, n q , see Fig. 4, we find that the global param- 
eter dependence can be described by the scaling relation 

Ak c (e) = Ae a , e = en s V% ( 2 ) 

The data fit gives A — 3.0, a = 0.64 for /j, = and 
A = 4.8, a = 0.68 for /i = e. This result can be un- 
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FIG. 4: Dependence of the shift of the critical point Ak c (e) = 

k c — fc c (e) on rescaled imperfection strength e = en g ^/n^ for 

e = 2 x 10" 5 (diamonds), 4 x 10" 5 (triangles) and 8 x 10~ 5 

(squares); open/full symbols are for /i = 0, 8 < n q < 13 and 

fi = e, 8 < n q < 11 respectively; fc c = 1.8. The dashed lines 

show the scaling relation 



derstood from the following arguments. According to 
[IH llflj the time scale tf , on which the fidelity of quan- 
tum computation is close to unity, is determined by the 
parameter e (tf ~ 1/?)- Thus, an effective matrix ele- 
ment induced by static imperfections between ideal lo- 
calized eigenstates can be estimated as U e f ~ (Q ~ e/Z , 
where Q is a typical overlap of localized eigenstates which 
for the Anderson localization in d dimensions can be es- 
timated as Q ~ Z _/3 with /3 = d/2 and Z being the local- 
ization length for the exact algorithm (see a discussion in 
|24| for cZ = 1). The imperfections induced delocalization 
should take place when U e f exceeds the level spacing in 
a block of size I (U e f > A; ~ 1/Z rf ). Taking into account 
that near the critical point the localization length scales 
as I ~ Ak~ v with v ~ 1.5 (see HlHi^) we obtain that 
a = \/(v(d — (3)) = 2/ ud. The obtained value of a would 
give a reasonable value of v m 1.0 but in our model 
the situation is more complicated. Indeed, the dynamics 
in Q takes place in one dimension and hence one ex- 
pects f} — 1/2 and v w 0.6. The later value has a notice- 
able difference from a usually expected value 0, 0, • 
A possible reason for this discrepancy can be related to 
the fact that in the algorithm the perturbations give far 
away transitions (see Fig. 1) which effectively decrease 
the value of /3, also near the critical point the correla- 
tions in the matrix elements can play an important role. 
Further studies are required to clarify this point. 

Finally, let us note that in the vicinity of critical point 
the number of states grows with time as n d ~ t 0, 0, Q • 
Hence, the number of classical operations for t kicks 
can be estimated as n ac i ~ tN d log d N ~ t 2 \og d t while 
the quantum algorithm will need n g ~ dr? q t ~ t log 2 t 
gates assuming d quantum registers with N d = 2 drlq ~ t 
states. The coarse-grained characteristics of the proba- 



bility distribution can be determined from few measure- 
ments of most significant qubits, e.g. W as in Fig. 3. 
Thus, even if each step in is efficient, the speedup is 
only quadratic near the critical point. Above the critical 
point we have diffusive growth with n d ~ t d / 2 and the 



(l+d/2) 



for d > 2. 



speedup is stronger: n gc i , tg 
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